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Abstract 

We obtain bounds for the spectrum and for the total width of the spectral gaps 
for Jacobi matrices on £ 2 (Z) of the form (Hip) n = a n _i^ n _i + b n ip n + a n ip n+ i, where 
a„ = a n+q and b n = b n+q are periodic sequences of real numbers. The results are based 
on a study of the quasimomentum k(z) corresponding to H. We consider k(z) as a 
conformal mapping in the complex plane. We obtain the trace identities which connect 
integrals of the Lyapunov exponent over the gaps with the normalised traces of powers 
of H. 



1 Introduction 

The purpose of the present work is to study the spectrum of the g-periodic Jacobi matrix 
on£ 2 (Z) 

{Hip) n = a n _i^ n _i + b n ip n + a n ip n+1 , (1.1) 

where b n = b n+q G K, a n = a n+q > 0, n G Z. The number q is the smallest period. (It is 

well known that any symmetric Jacobi matrix can be represented in the form with positive 

off-diagonal elements.) 

Periodic Jacobi matrices were discussed in many works (see, e.g., [H 121 El)- 

In what follows, we consider q > 1 thus excluding the trivial case of q = 1 (the spectrum 

of H with q = 1 is the interval [bi — 2a%, 61 + 20!]). Let r(H) denote half the distance between 

the ends of the spectrum of H. We set 

c = r(H), (1.2) 

A = (a 1 a 2 ---a q ) 1 l q . (1.3) 

We will also use the notation aj = a for a\ = a 2 = ■ ■ ■ = a q = a. We remark that periodic 
discrete ID Schrodinger operators are a particular case of (J1.1J) with aj = 1. 
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By Gershgorin's theorem (see, e.g., fIJ), the spectrum of H lies inside the interval 
[min,-(&j — dj — dj-i), maXj(6j + aj + a 3 -_i)]. It is absolutely continuous and consists of 
exactly q intervals 

0"m = [\n-i> \J, m = 1, . . . , q 

separated by gaps 

7m = (A~, A+), m = l,...,g-l. 

If a gap degenerates, i.e. 7 m = 0, then the corresponding segments a m , <J m+ i merge. In what 
follows, we will use | • | to denote the Lebesgue measure of sets. 

It is known [SJ |U| that the total width of bands and gaps, respectively, satisfy the in- 
equalities 1 

q AA q 9-1 

2 C -&>j> m |>— , £>mi>&, ( L4 ) 

m=l m=l 

where h = maxj bj — min^ bj and M = max(maxj(6j + aj + — mm,- bj, max,,- bj — minj(bj — 
dj — aj-i)). Analysis of the proof in [Hj shows that M in (|1.4)l can be replaced by 2c. 

In the present work we shall find further estimates for the widths of the gaps and for the 
total width of the spectrum in terms of the matrix elements of H. Our results are obtained 
from properties of the quasimomentum (associated with the periodic operator H) considered 
as a conformal mapping of complex domains. Originally, the method of conformal mappings 
was proposed in [Zj for the Hill operator (Hy(x) = y"(x) + f(x)y(x), f(x) = f(x + T)). It 
was further developed in [HI El UH1 HH H2j to obtain spectral estimates for the Hill operator 
and the Dirac operator. As we shall see below, the main peculiarity of the discrete case is 
that the spectrum is bounded. This makes z = arccos(A/ max|A^|) a more natural spectral 
variable than A. Thus our estimates for the gaps will be formulated in terms of z. 

As is known, the real part of the quasimomentum is the integrated density of states of 
the operator H while the imaginary part is the Lyapunov exponent (it describes the rate of 
growth in n of the solutions ip n (X) of the equation (H — \)tp = 0). In this paper we shall 
obtain identities (trace formulas and Dirichlet integrals) which relate various integrals of the 
quasimomentum to traces of powers of H ( Lemmas 13.21 and 13. 3|) . These identities will serve 
as a starting point to derive our estimates. They are also of separate interest. 

We shall call a Jacobi matrix normalised if — A^ = A^ > 0. Obviously, this can be 
achieved for an arbitrary matrix by adding a certain constant. Below we always assume that 
H is normalised. 

1 They were written in |H] for the Schrodinger case (dj = 1) but the arguments are easy to generalise. 
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We shall need the matrix L defined by 



bi d\ + id2 

d\ — ld 2 &2 



if 9 = 2, 



/ bi ai id q \ 
di b 2 d 2 

d 2 b 3 d 3 

S-2 b q _i d q _i 

\-id q a 9 _i b q ) 



if q > 2. 



(1.5) 

It follows from Floquet theory (see the next section) that c = r(-ff) > r(L). Hence, 
averaging over the eigenvalues of L gives the following estimate (lower bounds): 

c 2j >^TrL 2j , j = 1,2,... (1.6) 

Below we obtain a different type of bounds for c: 

Theorem 1.1 (Bounds for the width of the spectrum) Let H be normalised and q > 1. 
Then 

c>2A, (1.7) 
c 2 fl + ln -H- ^ > - Tr L 2 . (1.8) 

Remark: The equality sign in (jl.7|) happens if and only if q = 1. 

This theorem will be proved in Section 3. The first part of it, (jl . 7|) . will also be proved 
by a much simpler method which uses only general properties of polynomials fLemma I2.2j) . 
Note that the quantity A can be interpreted as the logarithmic capacity of the spectrum (see, 
e.g., (ED)- One might ask if the inequality (jl.7|) between half the diameter of a set and its 
logarithmic capacity holds in a more general context. For the Schrodinger case (dj = 1), the 
inequality (J1.7)) reduces to c > 2. It was recently shown in jHj that c > 2 for an arbitrary, 
not necessarily periodic, real nonconstant sequence bj, dj = 1. 

As will be clear from the next section, one can also derive further inequalities similar to 
(jl.8j) . They are analogous to (jl.6j) but better. We obtain from them lower bounds for c. 
Consider, for example, Harper's matrix (see JH] for a review): dj = 1, bj = 2cos(27ra;j + 6), 
a = p/q, p and q G Z, 9 G R. Obviously, c < 2 + max^ < 4. Simple calculation shows that 
TrL 2 = 4q. Substituting this value into inequality (jl.8|) and solving the latter numerically 
gives c > 2.41. Note that this bound does not depend on p and q. Therefore it follows from 
continuity that it also holds for irrational a. 

Let 

h+ = max -arccosh |y4~ 9 det(A — L)/2\. 

Aeu7„ q 

We obtain the following estimates for the gaps of H in terms of the matrix elements and c: 
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Theorem 1.2 (Bounds for the gaps) Let H be normalised and q > 1. Furthermore, let 
g n = (arccos(A~/c), arccos(A+/c)), where (A~,A+) is the n'th gap in the spectrum of H . 
Then at least one gap is open and 

^ . . \n{c/2A) \n{c/2A) n n . 

£ \9n\ > «^t-^ > «Tzkm> (L9) 



n=l 



\n{2c/A) 

Y [ cos 2 ^x>f ( Hc/2A) + 1 --^-Tr L A, (1.10) 
tiJ gn h + \ 2 4 2qc* J> 1 ' 

\n(c/2A) v^ t ,2 , , „x 



i.ir 



2c 



max 

where h + satisfies < h + < ln(^ + |^ S|=i < m 
Remarks: 

• Existence of an open gap implies that h + ^ 0. 

• Just like for (|1.8jl it is possible to derive further inequalities similar to (jl.lOj) improving 
the estimates. 

The plan of the paper is as follows. In the next section we recall some facts from Floquet 
theory and prove a lemma which guarantees existence of at least one open gap if q > 1 
and validity of (jl.7j) . Note that existence of gaps for the Schrodinger case also follows 
from (jl.4j) [6 . In Section 3 we construct the quasimomentum, obtain the trace formulas 
and prove Theorem 11.11 In Section 4 we show how our construction fits into the theory 
of the general quasimomentum developed in UHl El and recall some properties of the 
general quasimomentum. In Section 5 we use these properties together with estimates for 
the Dirichlet integrals (Lemma 13. 3|) to prove Theorem 11.21 and to establish further bounds 
on the maxima of the Lyapunov exponent ( Lemma 15. lj) . Note that the proof of the bounds 
fll.9j1 and (jl.K))) is simpler than that of (jl.llj) . To establish the former, we only need Lemmas 
OGH and EH 



2 Preliminaries 

Recall some facts from Floquet theory (see, e.g., [I]). 

For a fixed A we introduce two fundamental solutions <f)(\), 0(X) of the equation (%p = <p 
or if; = 6) (Hip(X)) n = \ip n (\), n = — oo, . . . , oo, by the initial conditions <f>o = 0, (pi = 1; 
#o = 1, #i = 0. Both of them satisfy the recurrence relation: 

a n if) n+1 (X) = (A - b n )if) n (\) - a n -i^ n -x{X), n = l,2,.... (2.1) 

The polynomial of degree q 

D(\) = ( j )q+1 (\)+e q (\) (2.2) 
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Figure 1: The discriminant is sketched for q = 3. The bands of the spectrum are shown by 
thick lines. 

is called the discriminant of H. 

For example, when q = 2 we have D(\) = (A 2 — {pi + b 2 )X + bib 2 — a\ — a|)/(ai0 2 )- 
It is known that: 

1) Any solution of the equation (H — X)ip = (except possibly for those corresponding to 
-D(A) = ±2) is a linear combination of two solutions ip + and ip~ with the property 

^± +? (A) = e ±iql ^il)±{\), fiGZ, where k(X) = iarccos^^. (2.3) 

q 2 

Henceforth, we shall fix the branch of arccos(x) by the condition arccos(O) = —n/2. The 
function k(X) plays a crucial role in the present paper. 

2) The spectral intervals o m are the image of [—2, 2] under the inverse of the transform -D(A) 
(see Figure 1). (Note that by ()2.3|) solutions ip n are bounded in n if -D(A) G [—2,2] and 
otherwise exponentially increase.) Hence, all the critical points of -D(A) are maxima and 
minima; moreover, at all maxima points A max -D(A max ) > 2 and at all minima points \mi n 
-D(Amin) — — 2- All q — 1 critical points are mutually different. 

3) The solutions 0(A) and #(A) satisfy the identity 

?+ i(A)0 ? (A) - <MA) e q+1 (X) = 1. (2.4) 
We shall need the representation of -D(A) given by 
Lemma 2.1 D(X) = A~ q det(A - L) for all A G C. 
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Proof. Denote Djj., j, k > 0, the determinant of the matrix A — L with the first j — 1 and 
the last q — k rows and columns removed, and L lq , L ql set to zero. Expanding det(A — L) in 
the elements of the first row, we get 

det(A-L) =D Xq - a\D 2q _ x . (2.5) 

Similarly, the expansion of by the last row gives: = (X—bj)Dij_ 1 —a?_ 1 Dij_ 2 , j = 
1,2, ... ,q, where -Di-i = 0, D w = 1. Comparing this recurrence with ()2.1|) . we get D\ q = 
a\a 2 ■ ■ ■ a q 4> q+ i by induction. Similarly, we show that D 2q -i = ° 1 °^" ag ~ 1 Q q . Substituting 
these expressions into ()2. 5|) proves the lemma. ■ 

Now we shall prove that q ^ 1 implies the inequality ()1.7|) and existence of gaps. The 
proof is based on the above given facts from Floquet theory and general properties of poly- 
nomials. 

Lemma 2.2 Let q > 1. Then c = r(H) > 2A and at least one gap is open. 

Proof. Suppose that c < 2A. Then by adding a constant to H, we first ensure that the 
spectrum of H (Spec(if)) lies within the interval / = [— 2A, 2A] = Spec(if )- Here H is 
the Jacobi matrix with a, = A and bj = 0. Denote by Dh(X) the discriminant of H. It is 
a polynomial of degree q with the coefficient A~ q of A 9 . Let us set all dj = A and bj = in 
it, and denote the resulting polynomial Dh {^)- It is the discriminant for Hq viewed as a 
matrix of period q. Since Spec(ifo) has no gaps, we conclude that maxA G / \Dh {\) \ = 2. 

Let A = — 2 A, X q = 2 A, and A^, j = 1, . . . , q — 1 be the critical points of D H (\). By 
our construction, they all lie in /. By Floquet theory, \Du(\j)\ > 2 and the signs in the 
sequence D H (\j), j = 0, . . . , q alternate. Thus 

D H (X q )>D Ho (X q ), D H (X q _ 1 ) < Djj^X^x), D H (X q „ 2 )>D Ho (X q - 2 ), ... (2.6) 

The polynomial S*(A) = D H (X) — D H() (X) is of degree less than q. However, as follows from 
(|2.6|) . it changes its sign on / at least q times. Therefore S(X) = 0. Thus the strict inequality 
c < 2 A is impossible, and c = 2 A happens only when Dh = Dh - AVe shall now demonstrate 
that the last identity implies H = Hq, i.e., the smallest period of H is 1. 

The method is essentially borrowed from [3]. Assume Dh = Dh (we shall omit the 
subscript). Hence A^ = Xj and m&x\ e i |-D(A)| = \D(Xj )| = 2. Let Uj, j — 1, . . . , q — 1 be 
the zeros of the (q — l)-degree polynomial 4> q (X). Since n (A), n = 1, 2, . . . are the orthogonal 
polynomials corresponding to H, all the Uj are simple and belong to the open set (— 2A, 2A). 
Using (JHI), we get \D(vj)\ = + q (i/j)\ > 2. Therefore v$ = Xf, \D(uj)\ = 2, and 

9 q {vj) = ±1, where the sign corresponds to that of D(vj), j = 1, . . . , q — 1. Since 6 q (X) is a 
polynomial of degree q — 2, it can be uniquely reconstructed from its values in q— 1 points Uj 
by the Lagrange interpolation formula. From the definition of the discriminant, we obtain 
<f) q+ i(X) = -D(A) — 6 q (X). On the other hand, we also know the monic polynomial (i.e., with 
the coefficient 1 of the highest degree) <p q (X) = Yi^Zii^ ~ note that it is D'(X) up to 
a factor. We shall now see that these two polynomials determine H. Note first that the 
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coefficient of the highest (n — 1) degree of <p n (X) is (aia 2 • ■ ■ 1 - Therefore we have the 

following recursion for the monic polynomials <f) n (X) = A n_1 + a™~ 2 A n ~ 2 + ■ • • + a°: 

n+1 (A) + (b n - A)0 n (A) + al_j> n - x {\) = 0. (2.7) 

All coefficients of this polynomial are zero. This gives 

b n = a n n ~ 2 - a n _! = ^r 3 -^; 2 -^^- 2 . (2.8) 

Assuming n+ i(A) and 4> n (X) are given, we then obtain the polynomial n _i(A) from (j2.7j) . 
Performing this procedure successively for n = g, g — 1, . . . , 1, we reconstruct the coefficients 
ax, ... , a g _i and b±, . . . , b q . Finally, a g is obtained from the highest-degree coefficient a of 
g +i(A): a q = (aa>i ■ ■ ■a (? _ 1 ) _1 . From the uniqueness of our reconstruction, it follows that 
H = H . 

Thus for q > 2 we have c > 2A. Suppose now that the spectrum of H has no gaps, 
i.e., |-D#(A)| < 2 on [— c, c] (assuming H is normalised). Then considering the inequalities 
between Dh(X) and Dh {X) now at the critical points of Dh (\) on [— c, c], we obtain as 
before S(\) = Dh(X) — Dh (X) = 0. This contradiction shows that there is at least one open 
gap in the spectrum of H. * 



3 Trace formulas and Dirichlet integrals 

Our first aim is to investigate the asymptotics of the function k(X) for large imaginary 
A. We shall see below^that it is convenient to introduce instead of A a new variable z 
such that the mapping k(X(z)) becomes asymptotically the identity. After determining the 
asymptotics of fc(A(z)), we construct the regions into which the upper half-plane is mapped 
by the functions z(X) and k(X). This information is then used to obtain the trace formulas 
and Dirichlet integrals for the quasimomentum. 

First, we set c = A~ = \X$\ (recall the normalisation of section 1). Then, in the variable 
( = A/c, the spectrum lies within the interval [—1, 1] and its boundaries coincide with — 1 
and 1. 

Using the expansion 



arccos x 



we get qk(X) = arccos(D(A)/2) = zlnD(A) + 0(1/A 2 ") as |A| -> oo. Here, by Lemma I2~T1 



lnD(A) = lndet(A-L) -In A" = Trln(A - L) - hiA q 

q ln A + Tr ~ = Q ln A ~ ^ ~X^' 

.7=1 J 
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Further, using 1)3.1)1 . we have for |£| > 1 



A c c v-^ / 2? 

!ffl- = dn — - + dn2( = dn — - + arccos ■ 



7 =1 ^ 



A * f-f V / 2j(2C) 2i ' 

Definition. The function 

= fc(— CCOS Z) + 7T, (3.2) 

where fc(A) is defined by ()2.3)1 . ()2.2)1 . is called quasimomentum. 
From the preceding discussion, we have: 

fc(^) = 2 + iQo + — 4—, z — ► zoo, (3.3) 



cos J 2; 
i=i 



where 



c: 



l/jc j ) i Tr L J if 1 < j < 2g and odd; 



g ° ln 2A' Qj y(l/j2P) (. -(l/jc^Tr^ if 2 < j < 2g and even. (3 ' 4) 
In particular, 

Q 1 = -TrL, g 2 = I--l^TrL 2 . (3.5) 

Investigating the mapping of the boundaries (as we explain below), we obtain the cor- 
respondence of the domains in A-, z-, and fc-planes under the conformal transformations 
z(X) = 7r + arccos A/c, k = k(z) as shown in Figures 2-4 (for q = 3). The notation for the 
real and imaginary parts of z and k is fixed as follows: z = x + iy, k = u + iv . The points 
A+, z+, fc+ correspond to each other under the mappings; so do the points A~, z~, k~. That 
the upper A-half-plane is mapped onto the half-strip by the function z(X) is easy to verify. 
We prove first the mapping of the boundary dS R onto OSr (for finite regions S R and Sr) 
and then let R — *■ 00. To obtain the mapping of the boundary of A domain onto that of k 
domain, it is convenient to represent k(X) + it in the form 

1 f DW dD(X) 



q J D{ x+) y/4-D(\)* 

and integrate (keeping in mind Figure 1) along the real part of dS R with vanishing semicircles 
above the ends of the gaps A „• . 

We now write these results more precisely. Let 

A+ = C + \U7 n , C+ = {A : ImA > 0}, 
Zl = C° + \Ug n , C° + = {z : < Rez < n,lmz > 0}, 
9n = (z~,z+) } n = l,...,q-l, 
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® 



% lq 



K /q 



K k;=k| k;=k + k"=7i 



Figure 4: fc-plane 

where z^ = tt + arccos A^ / c. The function z = tt + arccos A/c is a conformal mapping of A + 
onto Z+. The gap 7„ is mapped onto g n , n — 1, . . . , q — 1. Let 

/C° = C° \ Uc n , c n = (—, — + ih n ), n = l,...,g-l, 

where c n is a vertical slit. The function k = k(z) conformally maps onto KP + . The spectral 
interval [z^_ t , z~] is mapped onto the interval of length n/q of the w-axis [(n — l)n/q, nn/q], 
n = 1, . . . ,q; the gap onto the slit c n , n = 1, . . . , q — 1. The height of the n-th slit is 
h n = (l/g)arccosh |-D(A maXjn )/2|, where A maXjn is the critical point of -D(A) in the gap j n . 
Thus we proved 

Lemma 3.1 (quasimomentum) The function k(z) = k(-ccosz) + tt is a conformal map- 
ping of the domain Z\ onto the quasimomentum domain It possesses the asymptotic 
expansion \3. 

We shall now obtain the identities (trace formulas) which connect the integrals of u = Re k 
and v = Imk with the coefficients Qi in the asymptotic expansion of k(z) and, via these 
coefficients, with the matrix elements of H. The trace formulas will be the basis of all the 
spectral estimates in the present work. To establish (jl.7j) (in a way different from that of 
Lemma only the first of these formulas is sufficient. The other ones can also be used 
to extract some local information about the gaps, but we shall not do this here. 
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Lemma 3.2 (trace formulas) The following identities hold: 

i r c 

- v(x)dx = Q = \n—, (3.6) 

7T In LA 



1 



v(:c) cos n xcb = <j , n- \ ( 3 - 7 ) 

T • / " 1 y ) /2-* if n even. 



In particular, 



v(x)cosxdx = Qi, (3.8) 

1 r n 

— / v (x) cos 2 xdx = Qo/2 + Q 2 . (3.9) 

n Jo 

Proof. Let us calculate the integral of k(z) along the boundary of as shown in Figure 3. 
We integrate first along the contour OSr. 

p pir pR pO pO 

0= / k(z)dz = / k(x)dx + / (n + iv(n, y))idy + / k(x + iR)dx + / iv(0,y)idy. 
Jds R Jo Jo Jn Jr 

Now take the limit R —>• 00 using ()3.3|) for fc(a; + iR). We obtain 

/»7T pOO 

k(z)dz = / k(x)dx — iirQ — n 2 /2 — / (v(tt, y) — v(0, y))dy. 
Jo Jo 

Separating the real and imaginary parts, we obtain (|3.6|) and the identity 

r ti 2 r°° 

/ u(x)dx = — +l (v(n,y) - v(0,y))dy. (3.10) 
Jo 1 Jo 

Similarly, considering the integrals of (k(z) — z — iQ — iYTj=iQjl cos- 7 2) cos n 2;, n = 
1,2,... along dZ^ we get 

pTT n pTT 

/ v(x) cos n xdx = ~S^Qj I cos n_J xdx, (3-H) 
Jo j=0 Jo 

whence the lemma follows. ■ 
Proof of Theorem \l.l\ 

i) Since v(x) > 0, it follows from (|3.6|) that c > 2 A. 

ii) Similarly, the right-hand side of ()3.9|) is nonnegative. Substituting there (|3.4|) . we 
obtain (jl.8|) with the greater-or-equal sign. 

To make the inequalities strict, we use Lemma l2~2l if q > 1 there exists at least one open 
gap, hence v(x) is not identically zero. Thus the integrals in Lemma T3. 21 are strictly positive 
and the theorem is proven. ■ 
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Similarly, using ([3.7)1 . we can obtain further inequalities for c. 

For a function f(z),z G C + , we formally define the Dirichlet integral 



-ft \f{z)\ 2 dxdy. 



1(f) = ~ 

7T 

We shall now present two such integrals the first of which we use below for our estimates. 
Lemma 3.3 (Dirichlet integrals) The following identities hold: 

I(k(z)-z) = Qo = hi^; (3.12) 

I ([k(z)-z-iQ ] cos z) = 9^ + Q 2 -2Q Q 2 -^-. (3.13) 

Proof. Using the Cauchy conditions, we rewrite the integral I(f(z)) in the form 

/(/(*)) = - // \V\mf{z)\ 2 dxdy. (3.14) 

7T J Jz° 

Recall the Green formula for the harmonic function u>(x,y) = lmf(z) in the closed area Sr. 

\Vu(x,y)\ 2 dxdy = / uj—dl, (3.15) 



s R Jds R dn 

where d/dn denotes the normal derivative (the normal points outside Sr) to the curve 8Sr. 

Let us calculate ()3.15|) for / = k(z) — z, that is for u(x,y) = v(x,y) — y. We use the 
Cauchy conditions to see that the integrals along both vertical boundaries of Sr vanish. To 
calculate the integral along the interval [0, tt] of the real line, we note (cf. Figure 4) that 
u = const on each gap, and v = on each band, hence the term v(x)du(x)/dx vanishes 
on [0, 7r] and the contribution of this part of the path is L v(x)dx. Taking then the limit 
as R — ► oo (using the asymptotics of k(z) on the part of 8Sr where y = R to see that the 
contribution of that part to the integral tends to zero), we obtain 

\k'(z) — l\ 2 dxdy = / v{x)dx. (3.16) 
I Jo 

In view of (|3.6|) . we have (|3.12|) . 

We turn to the calculation of /(/) for / = [k(z) — z — iQ ] cos z. Proceeding as above 
and using ()3.9|) . we get 

I([k(z) - z - iQo] cosz) = -O + Q 2 + ^/', (3.17) 

4 Z7T 



P7T 

I'— (vu — vx + Q u) sin2xdx. (3.18) 
Jo 
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In order to evaluate I', we consider 



r ( k(z) 2 z 2 3 

lim / -V- " H z )* + iQ k(z) + — + -QI- t Q z+ 

2Q Qi , 2Q Q 2 + Ql/2 



cos z COS 2 z 



sm2zdz, (3.19) 



where d e Z^_ differs from dZ^_ only in the e— neighbourhood of the point z = 7r/2 where it is 
a semicircle above this point of radius e (the integrand has a simple pole at z = 7r/2). The 
integrand in ()3.19j) is chosen to satisfy two conditions: 1) integration of the imaginary part of 
the first three terms in it over [0, 7r] gives I' (j3.18jl : 2) the integrand behaves like 0(1/ cosz) 
as z — > zoo. We evaluate ()3.19|) as the contour integrals in the proof of Lemma 13.21 The 
imaginary part of ()3.19|) then yields 



Substituting this into ()3.17|) . we conclude the proof of the lemma. ■ 

Further formulas similar to ()3.121 l3~T3"j) can be obtained with more effort. 



4 Properties of the quasimomentum 

By Lemma f3-H the quasimomentum k(z) conformally maps the strip onto the strip K,° + . 
We can expand it to a conformal mapping of the upper half z-plane onto the upper half 
fc-plane (with slits) in the following way. First, let k(— x + iy) — k(x + iy) for x + iy e 
(reflection with respect to the y axis). Now continue k(z) periodically: k(z + 27m) = 
k(z), n = ±1,±2, .... Thus extended k(z) is a particular (periodic) case of the general 
quasimomentum: a well-studied mapping ^SJ We can utilise, therefore, the known 
properties of the general quasimomentum. The rest of this section will be devoted to a 
review of some of these properties. 
Define the so-called comb domain 

/c + = c + \r, r = ur n , r n = (u n ,u n + ih n ), h n ^o, nez, 

where u n is a strongly increasing sequence of real numbers such that u n — > ±oo as n — > ±oo, 
and {K}^ e £°°. 

A conformal mapping k(z) = u + iv from the upper half plane C + onto some comb /C + is 
called a general quasimomentum if k(0) = 0, and k(iy) = iy(l + o(l)) as y — >• oo. A general 
quasimomentum k(z) is a continuous function in z G C + . The inverse function z(k) maps 
each slit T n onto the interval (gap) g n = (z~, z+) of the real axis; and each interval u~\ 
onto (the band) [z^_ x , z~]. Obviously, v(z) > for z G C + . The function v(x) is continuous 
on the real axis, equals zero on the bands, is positive and reaches the maximum h n in the 
n's gap. Hence v G L°°(R). It is known that also u' x (z) = v' y (z) > for z G C+. On the real 
line the function u'(x) is positive on the bands and equals zero on the gaps. (Thus it plays 
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a similar role for the bands as v(x) for the gaps. However, it is more difficult to obtain good 
estimates for u'(x).) By the Herglotz theorem for positive harmonic functions (e.g., [T7]). 
we have 



V{z) = v ( 1 + Ul\f^ dt )' :£C+ ' (4 ' 1! 



For its harmonic conjugate: 



1 r°° ( i _ x i 
u ^ z )= x j r -j ( — _____ j v (t)dt + const , z G C+. (4.2) 
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(The term with ^h; is here only to ensure convergence if v(t) does not decay at infinity.) 
Further, 

1 f°° , x / 1 * \ , - 1 f°° v(t)dt , . 

where the value of the constant Co is obtained from the condition k(0) =0. 

Since by the definition of the quasimomentum v' y {iy) = 1 + o(l) as y — > oo, we have 
u' x {iy) = 1 + o(l) as y — > oo. Hence, the Herglotz representation for u' x has the form 

= 1 + - r ^r 2 dt, z G C+. (4.4) 

Consequently, on each gap g n , 

u'{t)dt 
J-oo (t - xf 

i.e. the function v(x) is concave on each gap. 

For our estimates we shall need the following result proved in |Hj: For any gap Qj one has 



- v 'L\ylo = (OyUo = ~ I ~ — > 0, x G g n , (4.5) 



V [X) = Wj{X 



+ 1 L raw)' aM = liz - zt)(z - z7)|1/2 ' x 6 *■ (4 - 6) 

Thus, the graph of the function v(x) lies above the semicircle over the gap (see Figure 5). 

Note that in our case the quasimomentum is periodic (with period 2ir). This allows us 
to reduce the general formulas. For example, (jUBj) takes on the form 

Lemma 4.1 Let k(z) = k(z + 2ir), z G C+. Then 

k(z) = Z + C + — / v(6) cotan I ) g?0, C = V.p. / v(6) cotan-d#. (4.7) 

27r y__ \ 2 J 2tt 2 

Proof. We represent ()4.3|) as a sum of integrals over the intervals [— 7r + 27m, 7r + 27m], 
n = — oo, ...,oo, change the integration variable in each of them to transform it to [— 7r,ir], 
and then use the periodicity of v(x) and the identity 

_ r ^ 1 1 

V.p. > = - cotan . 

^ ^ e-z + 2-nn 2 \ 2 

n=— oo x ' 

The value of the constant is obtained from the condition k(0) = 0. ■ 
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v(x) 




Figure 5: The graph of v(x) is sketched on a gap x G (z- , zf). 

5 Estimates 

We return to the quasimomentum we constructed in Section 3. Recall that for q > 1 the 
function v(x) is not identically zero. Recall also that we assume H to be normalised. 

First, as a side result, note that similarly to the case of the Hill equation jS], we have a 
double-sided estimate for Q in terms of the gap lengths \g n \ and maxima h n of v(x) in the 
gaps: 

Lemma 5.1 Let q > 1 and h n = max l£jn v (x). Then 

^ 9-1 ^ 9-1 

^-y^nbnl < Qo < -yZ h n\9n\- (5-1) 

n=l n=l 

Proof. We have 

1 r 1 9-1 

Qo = - v(x)dx < -}h n \g n \, 

which is the r.h.s. inequality of (|5.1|) . To establish the l.h.s. we observe that the concavity 
of v(x) on gaps implies that 

v(x)dx > —h n \g n \ 



and then use ()3.6|) . ■ 

In order to find a bound on the gap width from below, we shall need an upper bound on 



Lemma 5.2 Let q > 1 and h + = max n h n . Then 



I 4 • -}->> ) <hw. (5.2) 
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Proof. We begin with an argument similar to that in Let Xj be the eigenvalues of L 
(zeros of D(X)). By the inequality between the algebraic and geometric means 

q 

A"\D(\)\ = H\\-\ 3 \<(S(\)/qy, 

3=1 

where S(X) = Ylj=i \^ ~ Observe that for A G [min A n , max A„] 

S(X) < max < ^^(Aj — min A„), ^^(maxA„ — Xj) > . (5.3) 

I j=l 3=1 J 

Hence, 5(A) < cq + \ Ylj=i -M = C( i + I Sj=i fyl < 2cq, where we used the fact that H is 
normalised. Now note that 

h+ = -arccosh BMI < I \ n \D(X+)\ < \n(S(X+)/Aq), 
q I q 

where A + is a point of maximum of |-D(A)| in the gaps. The positivity of h + follows from 
Lemma 12.21 ■ 
We have 

i r i 9-1 h 9-1 

-/ v(x)dx < -)h n \g n \< —) \g n \, (5.4) 

71 I n 71 ^— ' 71 

1/0 n=l n=l 

and, in view of ()3.6|) . we have the lower bound ()1.9|) for the total gap width. The inequality 
()1.10|) is proved similarly using ()3.9|) . 

To prove the r.h.s. of we use the semicircle property of v(x): 

g = in- = -y «(x)dx>-x; — 2 — = 8^ KI ' 

1/0 n=l n=l 

which yields the r.h.s. of 

We are yet to prove the l.h.s. of For this we shall make use of the Dirichlet 

integral lETT^) . 



Lemma 5.3 We have the following estimate 

q-l 



f qh+ 1 

2_. h 2 n < 7r 2 6 + ln(c/2A), where b + = max < 1, — i > 

n=l ^ 71 ' 



(5.5) 



Proof. First, recall the following result from [TT]. Let a real function f(u + iv) belong to the 
Sobolev space W^D^h, a, j3)), where the domain D(h,a,ft) = {a < u < (3, < v < h} for 
some ft > a, h > 0. Let / obey the following conditions : 

(a) f{u + i0) = 0, if u e (a, ft), 

(b) / is continuous in the closure D(h,a,ft). 
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Then it was shown in ^1] that 



\f(+a + iv)\ 2 dv 7i f h \ it ,,..-> , , 

! 1 <-max<{l,- } II \Vf\ 2 dudv. (5.6) 

D(h,a,P) 



2 I /3 - a 



Now take the function /(w + iv) = lm(k — z(k)) = v — y(/c) and note that 

(a) the identity y(u + iO) = yields f(u + iO) = 0; 

(b) as the function y{k) is continuous in 1C\ , f is also continuous there. 
Therefore the conditions leading to (|5.6j) are fulfilled for this f(k) and the domain 
D(h n , 7in/q, 7i(n + l)/q). Since f(iin/q + iv) = v if < v < h n , we get 

= 2 J vdv ^ 7rmax Slhi^ J J \z'{k) — l\ 2 dudv. (5.7) 

Replacing max jl, — } by its maximum b + and summing over n, we obtain 

y^hl<7rb+[[ \z'(k) — l\ 2 dudv = nb + f f \k'(z) - l\ 2 dxdy, (5.8) 

where the last equality follows by change of variables. In view of (I3.12J1 we finished the proof. 
■ 

Applying Lemma IB~3l we get 

9-1 
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i /-7T i 9-1 j 9-1 9-1 I 9-1 

J ° m=l \ m=l \ m=l V \ m=l 



which yields the l.h.s. of (ll.llj) . 
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